Nonequilibrium quantum criticality in open electronic systems 
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A theory is presented of quantum criticality in open (coupled to reservoirs) itinerant electron mag- 
nets, with nonequilibrium drive provided by current flow across the system. Both departures from 
equilibrium at conventional (equilibrium) quantum critical points and the physics of phase transi- 
tions induced by the nonequilibrium drive are treated. Nonequilibrium-induced phase transitions 
are found to have the same leading critical behavior as conventional thermal phase transitions. 
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A central issue in condensed matter physics is the be- 
havior of systems as one tunes parameters (for example 
pressure, or magnetic field) so as to change the symme- 
tries characterizing the ground state 0, 0, H, 0| ■ The 
parameter values at which the ground state symmetries 
change (for example, from ferromagnetic metal to para- 
magnetic metal) define a quantum phase transition point 
(quantum critical point). At quantum phase transitions, 
spatial and temporal fluctuations are coupled, so that 
continuous quantum phase transitions in equilibrium sys- 
tems are typically described by critical theories involving 
an effective dimensionality d e ff greater than the spatial 
dimensionality d. 

While equilibrium quantum phase transitions have 
been extensively studied, the generalization to nonequi- 
librium conditions raises a largely open class of ques- 
tions. Nonlinear transport near a superconductor- 
insulator phase transition [|| and a ferromagnetic tran- 
sition driven by current flow in a closed one-dimensional 
system Q have been studied; however, a general system- 
atic understanding is lacking. 

In this paper we formulate a theory of nonequilibrium 
quantum criticality in itinerant electron systems coupled 
to reservoirs (c.f. upper panel of Fig. [1]) with which par- 
ticles may be exchanged. Nonequilibrium is imposed by 
differences between reservoirs; our systems are therefore 
subject to a time-independent drive, and are not charac- 
terized by any conserved quantities. A generic phase dia- 
gram is shown in the lower panel of Fig[TJ we take a sys- 
tem which at temperature T = may be tuned through 
an equilibrium quantum critical point by varying a pa- 
rameter S through a critical value S c and determine the 
changes induced by a nonequilibrium drive (generically 
denoted as V). Of particular interest is the transition 
generated by V if the V = system is ordered (vertical 
arrow in lower panel of Fig. [I]). 

Analysis of nonequilibrium systems proceeds from the 
time dependent density matrix p(t) defined in terms of 
a Hamiltonian H and an initial condition p(ti n u) via 
p{t) = e- lti{t - u ^ p[t imt )e iil ^- u ^ . The open systems 
we consider possess a well defined long-time state defined 
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FIG. 1: Upper panel: Schematic view of systems studied: in- 
teracting electron system coupled to two leads. Lower panel: 
Schematic phase diagram in a plane of equilibrium distance 
from criticality, 8, and departure from equilibrium V. The 
quantum critical point S c separates the long range ordered 
(8 < S c ) and disordered (8 > S c ) phases. The solid line de- 
notes a nonequilibrium phase transition. The dashed curves 
indicate the crossover from low V essentially equilibrium 
physics to the higher- V nonequilibrium-dominated regime. 



by a density matrix pss independent of the initial con- 
dition p{ti n it). We use the usual Suzuki- Trotter breakup 
to express the problem as a two-time contour functional 
integral [1, 0, [l(J, Hubbard-Stratonovich techniques to 
decouple the interaction terms by introducing auxiliary 
fields which we interpret as order parameter fluctuations, 
and integrate over the electronic degrees of freedom. The 
result is that the physics of the steady state system is ex- 
pressed in terms of a generating functional, Z, of source 
fields, 77,[1 E3 

Z(r)) = fV[mi,mf]pss,A(mi,mf)f V[m+(t),m-(t)] 

gS|fI{m+,m-,i)}] (1) 

Here m± are the fluctuating order parameter 
fields of interest and A is a short-distance cutoff. 
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/ T>[{m + (t),m,-(t)}] denotes an integral over all paths 
in function space beginning at nij on the + contour at 
t = and ending at mj at t = on the — contour. 
The contributions of paths with given endpoints are 
weighted by the steady-state reduced density matrix 
PA[{ m i(x),mf(x)}], whose diagonal elements describe 
the probability that at an instant of time the long 
wavelength components of the order parameter field take 
the configuration m(x). Eq [1] is the generalization to 
nonequilibrium systems of the partition function used to 
treat criticality in equilibrium systems. The important 
differences are the two time contours and the presence of 
pss which obeys a kinetic equation determined from 
the requirement that correlation functions calculated 
from EqQ]are causal and finite. Eq[T]may be analysed 
by the renormalization group method of integrating out 
modes near the cutoff and rescaling. 

We apply the formalism to the case of a two dimen- 
sional itinerant magnet placed between two noninteract- 
ing leads, with current flowing across the system (cf Fig 
[1]). We model the system by the Hamiltonian 



H Hlayer ~\~ H m i x -\- Hi eac [ s 

Hlayer ^ ^ id^.i-^g a ^i,a ~t~ Hint 



H,, 



(2) 
(3) 

(4) 



,k,a,b—L,R 



where the combination of the band structure implied by 
t and the interactions Hi nt is such that the isolated layer 
is ferromagnetic. The lead electrons are described by 

operators a and have free fermion correlators (ala^ — 
fb with a fermi function with lead-dependent chemical 
potential /i . For V b k ^ neither the number of electrons 
nor the magnetization in the layer are conserved. 

We assume, as is usual in studies of quantum critical 
phenomena [3, 0], that the electronic propagators and 
susceptibilities in the layer take the usual fermi liquid 
form and treat the interactions by a perturbative renor- 
malization group. The presence of the leads implies that 
the Greens functions describing the propagation of elec- 
trons in the interacting layer are 



= (G J 



ReY. R {p,uj)Y + \ImT. R (p^)\ 2 



(5) 
(6) 



where p is a two dimensional momentum within the layer, 
Ira's 11 — J2 a ^ a w ith F a (p, lj) being the rate at which 
electrons escape from the active layer into the lead a and 

Y> K = -ZiYsa^i^R^aiP'^) (! ~ 2 /( w - Ma)) determines 
the distribution function imposed by coupling to reser- 
voirs. We shall focus on excitation energies less than T a 
and momenta less than Im'E R /vF where nonconservation 
due to escape into the leads is dominant. The analysis 



sketched above then leads to a generating function of the 
form of Eq[fl with S K = S' (2) + S' (4) + •■■ where 



S(2) = -ijdtj dt' J d d r J d d r' (m cl (t, r) m q (t, r)) 

, AN 



o be 1 ]' 



m q (t',r<) 



(J) 



Here m„ 



,m c i = — j - ^- 



L q — g > "'ci — 2 ' em P s i s denotes 

terms of higher than fourth order in m and the fourth 
order term will be presented and discussed below. 
The quadratic-level inverse propagators are: 



.(8) 



( q , n) = -2i j: ab coth (9) 



The key quantity is Xl K . In equilibrium systems at 
T = 0, H K (t) vanishes at long times (as a power 
law for itinerant-electron models); however at T ^ 
and (for all of the models we have studied) out 
of equilibrium H K (t — > oo) ^ 0. Mathematically, 
Hk acts as a mass for the quantum fluctuations. If 
T or V = \fi a — 7^ 0, quantum fluctuations 
are gapped and at long times the theory is classi- 
cal. The (x" 1 ) ( ' R ' j4) describe non-conserved (because 
of the leads) overdamped magnetization fluctuations. 

{lab)- 1 = (T a (p,n = o)T b (p,n = o)/r 3 (p,n = o))\ FS 

are fermi surface averaged decay rates, and 7 _1 = 
Y^ a b=L R 7a6 ■ ^he " mass " (distance from criticality) 5 
depends on the interaction, layer density of states and 
coupling to the leads. The overdamped dynamics im- 
plies that even at V = T = the momentum conju- 
gate to m(q,t) is logarithmically large, so that the dc 
fluctuations are essentially classical. The density matrix 
is then easily obtained from the generating function us- 
ing the techniques of [l(J. We find, up to corrections of 
0(V 2 /r 2 ) in the argument of the exponential, 



p[mi(k),mf(k)] ~ 5 mi ,m f exp 



2Re[ X - 1 (k)] R \m l (k)f 



ijU K (n = o) 



(io) 

After calculation, we find that the leading nonlinearity 



is 



=-ij {d{k}) £ 



Uim q m cl 



(11) 



Here the Ui are interaction functions which depend on the 
momenta and frequency {k} of all of the fields. The level 
broadening due to the leads means any space-dependence 
may be neglected. In an isolated system with Hamil- 
tonian dynamics one would have u\ = u% independent 
of frequency and ^2,4 = 0. The coupling to a reser- 
voir means that the interactions are retarded and that 
ui.i 0. The limit of u± t 3 as all momenta and frequen- 
cies tend to zero is real and positive, so to obtain the lead- 
ing long- wavelength, low energy behavior we may treat 
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Mi. 3 as constants. However, at T = V = 0, U2,4 — > as 
the external frequencies are set to zero so the frequency 
dependence must be considered. This is somewhat in- 
volved, but the case important for subsequent consider- 
ations is when the two quantum fields carry a frequency 
±f2; in this case at V — 0, u 2 — ^ u 2 fi coth iy, while 
at T = 0, V 0, u 2 — » tig | V|, with real and positive 

%11j 2 , %%l 2 • 

We now formulate a renormalization group treatment, 
following along the usual lines [H, 0, IH : we choose a b > 1 
and in Eq [7] integrate out those fluctuations with mo- 
menta between A and A/6, treating the interactions per- 
turbatively. One technical remark is needed: causality 
implies (m q m q ) correlator vanish identically at all times. 
In a theory with a frequency cutoff, care must be ex- 
ercised to ensure that the cutoff does not violate the 
causality requirements. We find it simpler to work with 
a theory with a momentum cutoff but no frequency cut- 
off, so that we eliminate all frequencies for each removed 
momentum mode. We next rescale q, u>, T, V in order to 
keep the cutoff, coefficients of q 2 and ifl, and arguments 
of U K invariant; thus g— »^-,f2— s-p^-, T, V— > 2^1. 
Observe that under this rescaling the density matrix pre- 
serves the form Eq[l0l The result is a change in the mass 
8 



dS 
d In b 



= 2S + S Ul (b)g 



(12) 



A similar renormalization coming from u 2 preserves the 
form of the term but changes the coefficient. We 
interpret this as a finite renormalization of the broad- 
ening parameter 7 and do not consider it further. The 
mode elimination leads also to a renormalization of the 
interactions (note: U2,4 = *^2,4) 6 = 4 — d — z) 
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11*2,4, e 

■lttKR f fKA^ + l 2ui u 2 f RA (13) 



^- b =eu 2 - 2u 2 [15m (f KR + f KA ) - 2u 2 f RA ] 

+18ui [ Ul f KK - 2u 3 f RA ] (14) 
= c( Ul - u 3 ) ~ 12«! [6u 4 f RA - u 2 f KK ] + . (15) 

(16) 



^-=eu 4 -6u 2 u 3 (f« R + f« A ) + 



d la b 

Here g = «f J %xk; f ab 

R 



X°~ = X a (A, (J) is one of x > X 



= > / ^X a X h where 

= x R x A i an d 



Kd = J p ?j <5(g — A) (in our discussions we will set A = 
1). Note that the initial values are u\ — u 3 = 0(T 2 , V 2 ), 
u 2 ,U4 = 0(T,V), f KK = 2f RA + 0(T 2 ,V 2 ). The func- 
tions /, g depend on 8(T, V). When 8, T, V -> 0, f ab and 
g tend to constant values of order unity. For T, V > 1 
but 5 < 1, / KK ~ F 2 ,T 2 , (/ Kfl + / KA ),5 ~ V,T and 
~ 1. When 8 ~ 1 scaling stops. 
We solve the scaling equations starting from the 
physically relevant initial conditions 8, T, V -C 1, and 
taking d = z = 2 as appropriate for a thin mag- 
netic layer. To leading nontrivial order in T, V we 
need to retain only the first terms after the e term in 



= 21nb 



S + 3g J Q hlb dxui(e x )e 2x 

1+ Anoint > /l = 18 (f KR + f KA ) and the Sub " 

script denotes initial conditions. We integrate up to 
scales where 8 <~ 1 (initial conditions corresponding to 
region below dashed line in Fig[lJ or, if 8 remains small, 
max(V,T) ~ 1. For the latter case, we denote the value 
of 8 at the crossover scale by r. We focus on the inter- 
esting regime -l<r<0 corresponding to criticality or 
to a T = ordered state very near to the critical point 
(scaling trajectory depicted by vertical arrow in Fig[T]). 
We henceforth set T = 0. 

At the crossover scale the interactions are all small by 
powers of logarithms. In the crossover region V ~ 1 
the expressions are complicated. In the classical re- 
gion V ^> 1, the functions /, g acquire the V depen- 
dence noted above. Rewriting the scaling equations 
in terms of g = gVJ KK = f KK V 2 ,f KR + f KA = 
V{f KR + f KA )J RA = f RA ,vi = Ul V, v 2 = fia, v 3 = f 



Eqns [T3l [TU [T5l [TH We integrate these equations to 
obtain 8(b) 

«i(6) 



where 



and w 4 = ^ 



V'2 



leads to 



= 25- 



^5 = 2«! - 18x; 2 (/^« + /^) + 12« lVa / 
-^a- - 1 s^. 2 f KK - 30v lV2 {f KR + f KA ) 



RA 



-Avlf RA - 36 Vl v 3 f RA 



dV3 

din 6 



-2v 3 



dV4 

d\nb 



-4:VA 



(17) 
(18) 

(19) 
(20) 



Thus in the classical regime v 3 and U4 vanish rapidly, v± 
grows and v 2 reaches a fixed point. The effective the- 
ory becomes quadratic in m q which may be integrated 
out [l0(. The result is a theory for the fluctuations of 
the classical component of the magnetization field in the 
presence of a Gaussian, delta-correlated noise determined 
by the nonequilibrium drive: 

1 dm r i , q i-j „ s 

= (8-e o y 2 + vim 2 cl )m cl +^ 



7 dt 

with T eS 



= 8{x-x')8{t-t') 



vim cl )m c i 
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(21) 



1LR 



we have used a standard transforma- 
tion [llj to eliminate a coupling, generated by v 2l be- 
tween the noise and the classical field). All parameters 
acquire finite renormalizations (not explicitly denoted) 
from the scaling process. Eq[5l]is identical to that used 
in the standard analysis of Model A [l^] relaxational dy- 
namics, except that T e ff appears instead of temperature 
in the noise correlator. We therefore conclude that the 
voltage-driven transition is in the same universality class 
as the usual thermal 2d Ising transition, and more gener- 
ally that in this model, as far as universal quantities are 
concerned voltage acts as a temperature. From Eg 1211 we 
may extract experimental consequences: near the critical 
point, magnetic correlation length £ is = 8 ~ \ n X/v 
and the in-plane 2d resistivity arising from the scattering 
of electron with critical fluctuations is p(V) ~ V 3 ^ 2 . 
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We next extend the analysis to the 0(3) symmet- 
ric (Heisenberg) case. As in the equilibrium situation 
P, Q the physics of the disordered and quantum-classical 
crossover regimes is only weakly dependent on spin sym- 
metry. Differences appear in the "renormalized classical" 
regime corresponding to adding a weak noncquilibrium 
drive to an ordered state. In this regime the procedure 
leading to Eq. [21] gives a nonlinear stochastic equation 
describing fluctuations of the magnetization amplitude 
and precession of its direction. Here we focus on the 
most important special case, namely precession of small 
amplitude, low frequency, long-wavelength fluctuations 
of the magnetization direction about a state assumed to 
possess long ranged order directed along z and we denote 
the spin-gap by A. The important degrees of freedom are 
those transverse to the ordering direction. We find that 
at scales t > 1/7 and L > {vf/j) these are described by 



V T 



+ z x 



M; _ (h 
dt \ Ut 



byyAV 

77-L-R 



ZX £, 



(22) 



£ in Eq. [H] is a fluctuating noise field whose compo- 
nents are independent and correlated according to Eq 
1211 but with 7 replaced by 7/a X x. The a,b coefficients 
are numbers of order unity. The subscripts indicate 
whether they arise from the xx or xy terms in the re- 
tarded/advanced susceptibilities. The term involving a xy 
gives the Landau-Lifshitz-Gilbert spin precession and is 
explicitly proportional to the spin gap A so is small near 
the quantum critical point. The term involving a xx ex- 
presses the damping due to coupling to the leads, remains 
nonvanishing as A — ► and is the dominant time deriva- 
tive term. Solving Eq[H]in the rotating frame leads to 



(— ^- <C 1) there is no instability. Calculation reveals 
that obtaining a nonzero b xy also requires an energy- 
dependent asymmetry between the leads (Tl(si)^r(£2) — 

We have presented a theory for nonequilibrium phase 
transitions in an itinerant-electron system coupled to ex- 
ternal reservoirs. We provide a precise mapping onto 
an effective classical theory which demonstrates that the 
leading effect of the nonequilibrium drive is to gener- 
ate an effective temperature and hence a transition in 
the standard Wilson-Fisher thermal universality class. 
Nonequilibrium- induced breaking of time reversal and in- 
version symmetries and the creation of a coherently pro- 
cessing ( "spin-torque" ) state appear only at the level of 
subleading corrections. The techniques introduced here 
can be applied to important open problems such as sys- 
tems where the drive couples linearly to the order param- 
eter [jj Q , driven Bose condensates [l5[ , and the closed 
system, conserved order parameter work of 0]. 

Acknowledgements: This work was supported by NSF- 
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2005-070-C00044 (AM,ST,YBK) 

Note Added: Shortly after our manuscript was submit- 
ted a study of the closely related problem of nonlinear 
transport at a quantum critical point appeared [l6j |. 



(m+(q, t)m-{q't')) 



Eq. 



shows that the mean square magnetization 



fluctuations diverge as -3. This divergence signals the 
instability of the ordered state by the voltage-induced 
decoherence in precise analogy to the usual 2d ther- 
mal case. The oscillatory term in Eq [53] expresses the 
usual spin precession with precession frequency u; e ff = 



— )A<7 2 £q shifted from the equilib- 



rium result by an amount proportional to V due to spin 
accumulation effects at the interface of the magnetized 
layer and leads [T3| • The decaying term in Eq. [23] ex- 
presses the damping due to coupling to leads D c g — 



7) 



If 



< 



a b Av 

U xy U x y lLRl 2 



Eq [22] supports modes which grow exponentially with 
time leading to the spin-torque instability recently dis- 
cussed [14]. However, where the present theory applies 



I g( <? + g , ) e -°°H < ' 2 ^''- t ''-'^ff(*-;;u 
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